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Abstract 

We discuss the time development of Gaussian wave packet solutions of the 
'quantum bouncer' (a quantum mechanical particle subject to a uniform down- 
ward force, above an impermeable flat surface). We focus on the evaluation and 
visualization of the expectation values and uncertainties of position and momen- 
tum variables during a single quasi-classical period as well as during the long 
term collapsed phase and several revivals. This approach complements exist- 
ing analytic and numerical analyses of this system, as well as being useful for 
comparison with similar results for the harmonic oscillator and infinite well cases. 
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I. Introduction 



With the advent of modern computer technology, robust numerical calculations of 
time-dependent phenomena in quantum mechanics are now common, as is the visual- 
ization of the resulting effects [1] - [4]. The time-evolution of wave packet solutions for 
many scattering geometries as well as for bound state systems [5] are discussed with 
increasing frequency in the pedagogical literature, illustrating not only such familiar 
aspects as wave packet spreading, but also extending student experience to more novel 
phenomena such as wave packet revivals. Reviews of wave packet revivals (initially 
highly localized quantum wave packet solutions which exhibit quasi-classical behavior, 
then disperse or spread in time to a so-called collapsed phase, only to reform later to 
something very much like it's initial state) in many familiar model one-dimensional 
quantum mechanical systems such as the harmonic oscillator and infinite well have ap- 
peared [6] - [9] providing students with accessible examples of an important quantum 
effect which can be probed experimentally in atomic systems. 

Besides studies of such effects in the two 'classic' quantum mechanical model sys- 
tems mentioned above, wave packet propagation, and especially the structure of re- 
vivals, has also recently been discussed in the quantum version of another familiar 
classical system, the so-called quantum bouncer [10], [11]. This system is the bound 
state version of the classical 'falling object', and is defined by the potential energy 



corresponding to a constant downwards force acting on a particle above an impene- 
trable flat surface at 2; = 0. Generalizing on many early papers [12] which discuss the 
time-independent solutions for this problem, the author of Ref. [10] focuses attention 
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for 2; < 
for < 2; 



(1) 



on deriving expressions for the collapse and revival times for initially Gaussian wave 
packet solutions, using a mixture of numerical and analytic techniques. In this note, 
we will revisit this problem, along much the same lines as in Ref. [9], by examining, 
in detail, the short-term (quasi-classical) and long-term (revival) time-development of 
such solutions in terms of their position and momentum expectation values and un- 
certainties. This type of expectation value analysis, coupled with existing analytic, 
numerical, and visualization studies can then help form a more complete picture of 
the highly non-trivial time-development possible in one of the 'classic' one-dimensional 
model systems in introductory quantum mechanics. 

As often happens, the more familiar and tractable cases of the harmonic oscillator 
and the infinite well examples exhibit special features, in this case, in terms of their 
time-dependence, compared to a more general form such as in Eqn. 1. For that reason, 
in Sec. II, we first very briefiy review Gaussian wave packet solutions in these two more 
familiar systems. The bouncer system shares with the infinite well the feature of the 
'bounce' at the impenetrable wall, but some observables, such as the time-dependent 
wave packet spread, exhibit a much more typical cyclic structure than for the infinite 
well, where only the free-particle spreading time [9] is of relevance. Since the numerical 
studies presented here can also make contact with analytic solutions to the problem of 
a particle undergoing uniform acceleration, we also mention those solutions in Sec. Ill 
and then proceed to detailed results for the quantum bouncer. 

II. Gaussian wave packet solutions for the harmonic oscillator 
and infinite well 

Because of the special nature of the energy level structure of the harmonic oscillator, 
any time-dependent solution of this problem, ip{x, t), will be exactly periodic and return 
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precisely to its initial state after the classical period, Td = 27i/uj. Using standard 
propagator techniques or other methods, it is, in fact, easy to construct closed form 
Gaussian wave packet solutions [13] such as 

\ilj{x,t)\^ = -lx-xocos(u.t)]yLHt) 

\(p{p,t)\^ = _^__e-[P+~«i°('^*)lVpi(*) (3) 

which have expectation values which satisfy the classical equations of motion, namely 

{x)t = XQCOs{u!t) and {p)t = —mu!Xosm{u!t). (4) 

The fact that the expectation values for any wave packet solutions for the harmonic 
oscillator satisfy the classical equations of motion does not depend on these specific 
Gaussian forms, but can be shown to be true analytically in a quite general way 
[14]. The solutions in Eqns. (2) and (3) have time-dependent position and momentum 
uncertainties given by Axt = L(t)/\/2 and Ap^ = px,(t)/\/2 where 

[L{t)f = L^cos^{ujt)+(^^\ sin^M (5) 

\pL{t)f = cos^ {ujt) + {Lmujfsm^{ujt) (6) 

and the parameter L sets the scale for both the position- and momentum-space spreads. 



Unless L — ^Jh/muj, these position- and momentum-uncertainties oscillate with a 
period twice that of the classical motion [13], regaining the initial values of Axq and 
Apo at two opposing points in phase space. Once again, these results arc not specific 
to the Gaussian solution, but have been derived quite generally in Ref. [14]: they are 
also sometimes rediscovered [15], [16] in these specific cases. 

Another tool which is standardly used to probe the wave packet's approximate 
return to the initial state is the auto-correlation function [17], defined by 

/+00 r+oo °Q 

r{z, t) i;{z, 0)dz^ <f>*{p, t) <t>{p, 0)dp^Yl \^n\' e^^"*/^ (7) 
-oo J—oo n=l 



which measures the overlap of the position- or momentum-space wavefunction at later 
times with the initial state. For the oscillator, one can argue that because the wave 
packets are exactly periodic and never collapse, there are no revivals. 

On the other hand, for the infinite well, because the energy levels are integral 
multiples of a common value (but not equally spaced) , there are exact revivals [6] , but 
initial Gaussian-type wave packets do undergo dispersion into a collapsed phase. For 
example, wave packets constructed from energy eigenstates via 

oo 

^(a;,t) = $:c„«„(a;)e-^"*/'^ (8) 

n=0 

where Un{x) = sm(mvx / L) / \/L for the infinite well, with quasi-Gaussian expansion 
coefficients of the form 



(where p„ = nhir/L) give an initial momentum distribution very close to the standard 
one for Gaussian free-particle packets, namely 

and yield localized Gaussian-like position wave packets which are initially very close 
to the free-particle form 

\ij{x,t)\^ = ^_e-(^-N+P0t/m])V/3f (11) 



where Pt = ha^l + t^/tQ with to = inha^. Such packets undergo quasi-classical 
motion, bouncing back and forth between the two infinite walls, with an increasing 
width given roughly by an envelope defined by the time-dependent free-particle width, 
Axt = Pt/V^ [9]- During the impulsive collisions with the walls, the position-space 



width temporarily decreases during the time that the momentum-space distribution is 
'flipping' from positive to negative values. 

As mentioned above, since the quantum bouncer has features in common with both 
systems, namely one inflnite wall boundary at which 'bounces' will occur and a second, 
smoother potential barrier, we will keep both these cases in mind as we examine the 
short- and long-term behavior of expectation values for the quantum bouncer. 

III. Expectation values for Gaussian wave packets for the 
quantum bouncer 

Before turning to wave packet solutions of the quantum bouncer problem, we first 
recall some results from the related problem of a quantum particle undergoing uni- 
form acceleration, namely subject to a constant force, F, with a potential given by 
V{x) = —Fx everjrwhere in space. Transforming the resulting Schrodinger equation 
into momentum-space, one can construct arbitrary solutions [18] of the form 



where 0o(p) is any initial momentum distribution. This result already implies that 
the momentum distribution for this case simply translates uniformly in time, with no 
change in shape, since 



so that Apt = Apo- (We expect to see this behavior initially in the quantum bouncer 
case, until the wave packet nears the wall, as well as after the collision, at least for the 
first few classical periods.) Using this momentum-space solution, we can also construct 
closed-form position-space solutions; for example using the distribution 





(13) 




(14) 
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we find a general Gaussian wave packet for tfie accelerating particle 

\^(x,t)\' = A^e-^-\-^^V.tlm^FeilM?ie!l (15) 



with /3( = ah\Jl + t^/to- (The same results can be obtained in a variety of ways, using 
operator [19] or propagator [20] methods or other techniques [21].) This solution gives 

{x),^x^ + ^-r^e and Aa;, = ^A = ^Vl + ^'Ai- (16) 

m Zm V2 V 2 

Thus, the expectation value tracks the classical trajectory, while the wave packet 
spreads in exactly the same way as the standard free-particle Gaussian and we ex- 
pect to see remnants of this behavior as well. 

The author of Ref. [10] has constructed quasi- Gaussian wave packet solutions for 
the quantum bouncer using a combination of numerical and semi-analytic techniques, 
focusing on the derivation of quantities such as approximate revival time, Trev (See 
Ref. [6] for a very general derivation of both the classical period and quantum revival 
time for arbitrary bound state energy spectra). For an initial wave packet which is 
constructed so as to represent a particle released from rest at a height zq, the classical 
period of motion, Td, and the revival time, Trev, are given by 



r. = V^ and T„ = i[?^) (17) 
respectively. For simplicity, following Ref. [10], we take units such that 

h^l, 2m = 1, and F=l (18) 

and we choose particular initial values, namely 

zo = 25 and A^o = 1- (19) 
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With these values, the classical period and revival time are given by 

4 

Tci = 10 and Trev = -(5)^ = 795.78. (20) 

TT 

The collapse time, the time it takes for the initially localized wave packet to de-phase 
and become spread over the entire well, is given, in these dimensionless units as well 
as for our specific wavepacket parameters, as 

Tcoll = or Tcoll ^125 (21) 

respectively. 

Using these parameters, we can evaluate both the position-space and momentum- 
space wavefunctions in a manner which is similar to that in Ref. [10]. We generate the 
expansion coefficients, c„, in Eqn. (8) by explicit numerical calculation of the overlap 
integrals of the initial Gaussian form with the normalized energy eigenstates obtained 
by numerical solution of the Schrodinger equation; the eigenstates, normalization con- 
stants, and expansion coefficients are then compared to the Airy function solutions 
used in Ref. [10], supplemented by the more exact analytic results in Ref. [11], as 
a cross-check. (We truncate the expansion in eigenstates at a level consistent with 
the double-precision accuracy of our numerical evaluation, typically at values of the 
Cn ~ 10~^; this is to be be compared to the maximum values of c„ ~ 0.5. The resulting 
wavefunctions are then found to be appropriately normalized to unity to within 10~^ 
or better over all time intervals considered here.) 

We use this more numerical procedure partly because this technique runs more 
efficiently in our computing environment, but also because we intend to eventually 
extend our investigations of wave packet behavior to other, more general power law 
potentials of the form V(fc)(a;) = Volx/al*^ for comparison to recent work on the revival 
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times in such potentials [22]. This class of potentials contains the harmonic oscillator 
{k = 2) and the infinite well {k = oo) as special cases and can also be analyzed with 
an infinite wall added at x = to give the quantum bouncer case (where k — 1): this 
fact also helps motivate our brief review of those two special cases. 

We then plot the corresponding probability distributions over the first classical 
period in Fig. 1. On the left, the position-space probability is seen to spread in a manner 
which is numerically consistent with Eqn. (15), while the calculated position value {z)t 
(solid curve) agrees well with the classical expectation for the trajectory (dashed) 
except, of course, for the cusp at the 'bounce'. The packet exhibits the standard 
'interference' pattern during the collision with the wall [5], [23], at the 'bounce', and 
then reforms into something like the initial packet (compare to the dotted initial packet 
superimposed on the t = 10 case), only wider. 

For the momentum-space distributions (shown on the right of Fig. 1), we also see 
features of both the classical motion and the uniformly accelerated wave packet. The 
expectation value of momentum {p)t, calculated from |0(p, t)p and plotted as the solid 
curve, is once again consistent with the classical trajectory (dashed curve), except 
near the discontinuous, impulsive change in momentum values at the 'bounce'. The 
shape of the momentum-distribution follows the form expected from Eqn. (13), namely 
uniform translation with no change in shape, from t = 1 — > t = 3 and then again from 
t — 8 ^ t — 10, that is, during the time when it is not in collision with the wall, but 
with a definite final change in shape, compared to the initial |0o(p)P superimposed 
on the t = 10 result. The dotted vertical lines indicate the values of p = 0, and 
also the classically expected minimum and maximum values of momentum given by 
±Pm = ±V2mE = ±V2mFzo. 

The asymmetric shape of |0(p, t)p which is obvious after one period can be un- 
derstood using purely classical arguments. For example, we can examine the classical 
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trajectories of particles released from the same height, i7, but with varying initial ve- 
locities (or momenta); the resulting time-dependent momenta, p{t) versus t, are shown 
in Fig. 2. We can see the same asymmetry arising at t — Td ets for the quantum 
mechanical [^(p, t)P in Fig. 1; trajectories with momentum values p at t — are 
shifted upwards towards p ^ at t — T^, while those with p > are also moved to 
higher values to form the observed high-momentum 'tail'. (A very similar classical ex- 
planation of the time-dependent behavior of the momentum-space probability density 
for an otherwise free particle undergoing a 'bounce' from an infinite wall is discussed 
in Ref. [5]; in that case, the observed asymmetry in \(f){p, t)\'^ during the 'bounce' is due 
to the fact that the fastest (highest momentum) components of the wavepacket strike 
the wall first, and hence are reversed in direction, before the slow components.) 

We next focus on the time-development of quantum mechanical expectation values 
over the first eight classical periods and plot {z)t and {p)t over this interval in Fig. 3. 
Once again, the solid curves are the quantum mechanical expectation values, while 
the dashed curves are the classical trajectories. Clearly, the quantum results generally 
track the classical values, but become increasingly smoothed out as the packet spreads 
in time. We show in the same figure the time-development of the autocorrelation 
function, |^(i)P versus t, over the same time interval which shows the (increasingly) 
approximate nature of the return to the initial state at multiples of the classical period: 
these results look very similar to data for wave packets in the infinite well [6], [9] as it 
approaches the collapsed state, while the same plot for the oscillator gives a completely 
periodic structure with A{t) returning exactly to 1 at integral multiples of r = 27r/a;. 

In Fig. 4, we plot the time-dependent uncertainties in the position {Azt, top) and 
momentum (Apt, bottom) over the first eight classical periods. We note that the 
behavior of Az^ is somewhat similar to that for the oscillator case: it initially increases 



9 



in a manner consistent with Eqn. (16) for much of the first half period, but then shrinks 
dramatically near Tci/2 during the collision with the wall (as in Ref. [5]), and finally 
shrinks more slowly during the second half period, returning to something close to its 
initial spread. This behavior is different than the case of the infinite well where the 
wavepacket spreading is mostly determined by the free-particle expansion and more 
similar to the harmonic oscillator case in Eqn. (6) and general bound state systems. 

In Fig. 4, the expectation value during each period is superimposed with two ex- 
pressions for the spread (shown as pairs of dashed curves); specifically, for the first 
period we use 



and we extend these cyclically to later periods as shown in the figure. These agree 
fairly well with the observed timc-dcpcndence during the first cycle except during the 
collision phase, but become increasing bad representations as the wave packet spreads 
into its collapsed phase. Thus, over at least the first few periods, the return to the 
initial spread is similar to that in Eqn. (6). 

In the same figure, the spread in momentum is also illustrative of both classical 
motion and results for accelerated wave packets. Over the first cycle, when the packet 
is not colliding with the wall, the spread in momentum remains constant, as expected 
from Eqn. (13), while it increases dramatically during the collision time, just as with 
the 'bouncing wavepacket' in Ref. [5]. After each bounce, however, the spread in 
momentum has grown slightly so that Apt increases in a quasi-stepwise manner between 
collision, with flat plateaus between each impulsive 'spike'. 

We now turn our attention to the long-term behavior of these solutions and plot in 
Fig. 5 the same quantities, {z)ti {p)t, and |A(t)p, as in Fig. 3, but now over an interval 




and 




(22) 
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containing two expected revivals; these are indicated by the bold vertical dashed lines. 
We also indicate the collapse time, given by Eqn. (21), after t = as well as one unit 
of collapse time around each expected revival. The reformation of the wavepacket at 
the revivals in clearly apparent in all three variables, as is the approach to the more 
constant values during the collapsed phases. 

The horizontal dashed lines around which the quantum mechanical expectation 
values cluster during the collapsed phases are given by the average values derived from 
purely classical probability densities. Using standard arguments about how much time 
a particle spends in a small interval of position, we can find the normalized position- 
space classical probability as 



and which vanishes for all other values of z. The classical turning point, /I, is given 
by £^ = FA, and we naturally associate A with the initial position zq in our quantum 
analysis. With this identification, the classical average values are 



For the momentum-space probability distribution, we note that for a constant force, 
the probability density is uniform over the allowed space of values, namely 




for < 2; < A 



(23) 




(24) 



Pcl{v) 



1 



for 



(25) 



2pM 



-PM <P< +PM 



11 



and vanishing for \p\ > pm- in this expression we have pm = V'^mE = y/2mFA — 
\/2mFzQ. The classical averages are then given by 

(p) = 

(/) = (26) 



Ap = ^ (p^) - (p)^ = :^pm- 

These results are known to agree well (after local averaging) with the position- and 
momentum-space probability densities corresponding to time-independent energy eigen- 
states for this type of potential [24] . We would expect that the quantum wave packet 
expectation values to closely match these classical results (again, after local averaging) 
during the collapsed phase as well, since the wave function is closer to being an incoher- 
ent sum of many such stationary states rather than the highly coherent, well-localized, 
quasi-periodic state close to the initial time and later revivals. The time-dependent 
values of {z)t and {p)t in Fig. 5 do collapse to just these values. 

For the autocorrelation function we have no classical analog, but we can formally 
evaluate |^(i)P in terms of the expansion coefficients, c„, of the wave packet. In the 
limit that the various time-dependent components of the packet get out of phase, we 
have the expectation that during the collapsed interval 



H \Cn 
n=l 



2 ^iEnt/n 



E |cnnc^re^(^"-W^;^|ej4 (27) 

n,m=l n=l 



as the 'off-diagonal' terms tend to cancel as they have lose the phase coherence built 
into the initial state. This limiting value, calculated using the numerically determined 
Cn, is plotted in Fig. 5 (as the bold horizontal dashed line) and compares well to the 
locally averaged value of |>l(t)|^ during much of the collapsed period. 

In Fig. 6, wc plot the position- and momentum-spreads, as in Fig. 4, but once again 
over the longer time interval, and superimpose the classical (constant) expectation 
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values for A.z and Ap as bold horizontal dashed lines. In both plots, the quantum 
wave packet expectation values and uncertainties cluster appropriately around the 
classical results of Eqns. (24) and (26) during the collapsed phase as expected. 

Finally, one of the most interesting results of Ref. [10] is that the wave packet 
revivals happen in such a way that the probability densities are almost completely out 
of phase with the classical motion. In order to visualize this effect in the expectation 
value approach followed here, we plot again the same quantities as in Fig. 3, namely 
{z)ti {p)t, and |74(i)p, also over a time interval corresponding to 8 classical periods, but 
now centered around the expected revival time. In each case we indeed see that the 
quantum expectation values (solid) are approximately half a cycle out of phase with the 
corresponding classical trajectory values (dashed) and that the autocorrelation function 
is peaked at regular intervals, but offset from integral multiples of the classical period 
by roughly T,i/2. 

In conclusion, we have reexamined the interesting problem of quasi-Gaussian wave 
packets for the quantum bouncer, focusing on the numerical calculation and visual- 
ization of the results through the expectation values and uncertainties of position and 
momentum. We have focused on a single classical period to illustrate the important 
similarities and significant differences between this case and well-studied periodic wave 
packet solutions to other familiar problems such as the harmonic oscillator and infi- 
nite well as well as to the closed-form solutions for the constant acceleration case in 
Eqns. (15) and (16). We have also examined these quantities over the first few periods, 
well before the collapsed phase and compared them to classical trajectory expectations. 
Finally, we have confirmed, numerically and through visualization, that the quantum 
expectation values approach (in a locally averaged sense) the purely classical results 
from the time-independent probability distributions in Eqns. (24) and (26) during the 
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collapsed phase as well as illustrating the expected phase relationships between the 
classical trajectories and quantum expectation values near the revivals. 
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Figure Captions 



Fig. 1. Gaussian wave packet solutions for the quantum bouncer in position-space {{il^iz, ^ 
versus z, left) and momentum-space (|0(p, t)p versus p, right) for various times 
during the first classical period. The solid curves represent the time-dependent 
expectation values of position {{z)t, left) and momentum {(p)t, right) for these 
solutions. The similar dashed curves are the classical trajectories, z(t) (left) and 
p{t) (right), superimposed. The wave packet parameters in Eqn. (19) are used. 
For the momentum-space figure, the vertical dotted lines represent the values 
p = and the classical extremal values of momentum, ±Pm = ±\/2mFi^. 

Fig. 2. Plots of the classical momentum versus time, p{t) versus t, for objects released 
from a height H with initial momenta given by p(0) — (solid), p(0) = ±Ap 
(dotted), and p(0) = ±2Ap (dashed). After one classical period, the distribution 
of momentum values is asymmetric about p = 0, in just the same way as the 
|0(p, t)P plot in Fig. 1; compare the momentum values defined by the t = 0, 
{+2Ap, — 2Ap) 'band' to the corresponding values at t — T^. 

Fig. 3. The sohd curves represent the expectation values of position {{z)t, top) and mo- 
mentum {{p)t, middle) and the autocorrelation function (|74(t)p, bottom) versus 
time for the first eight classical periods of motion for the bouncing wave packet. 
The wave packet parameters in Eqn. (19) are used. The dashed hues are the clas- 
sical trajectories for comparison. The horizontal line in the top figure corresponds 
to the average value of position over one classical period, namely {z) = 2zo/3, 
from Eqn. (24). 
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Fig. 4. Same as Fig. 3 except that the time-dependent values of the uncertainty in posi- 
tion (A^t, top) and momentum (Apt, bottom) are shown over the first eight 
periods. The dashed curves in the top picture correspond to the spreading 
(anti-spreading) of a Gaussian wave packet as it 'shdes' downhill (uphill) as in 
Eqn. (22); these expressions agree well with the observed values of Azt during the 
first period, except during those times when the packet is 'bouncing' against the 
infinite wall at 2; = 0, near Tci/2 — b. The dashed line in the Ap^ figure (bottom) 
indicates the classical spread in momentum defined in Eqn. (26), Ap = VmI^P^- 

Fig. 5. Same as Fig. 3, but now over a longer time scale, containing two revivals. The 
expected locations of the revivals, according to Eqn. (17), are shown as vertical 
dashed lines, and times within one collapse time, Tcou^ of the initial time and these 
two revival times are also shown as vertical dotted lines. The classical average 
values of position and momentum ((z) = 22;o/3 and (p) = 0) from Eqns. (24) 
and (26) are shown as bold horizontal dashed lines, as is the expected value of 
the autocorrelation function during the collapsed phase, |^(t)P = Yl,n knl*^) from 
Eqn. (27), using the numerically determined values of the c„. 

Fig. 6. Same as Fig. 3, expect for the long term values of lS.Zt and Ap^. The classical 
average values of these quantities (A^ = 2zq/\/^ and Ap = Pm/V^) from 
Eqns. (24) and (26) are shown as bold horizontal dashed lines. 

Fig. 7. Same as Fig. 3, except for an eight classical period time window around the first 
revival which is indicated by the bold vertical dashed line. Note how the quantum 
mechanical averages, {z)t and {p)t (solid curves), are out of phase with the clas- 
sical trajectory values (dashed curves). In the bottom figure, the dashed curve 
corresponds to the time-development of the auto-correlation function, |A(i)p, 
over the first eight periods of the motion, for comparison. 
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